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Abstract 


A  methodology  is  presented  for  the  construction  of  process  models  by  the  combination 
of  physically  based  mechanistic  modeling  and  statistical  experimental  design  in  order  to 
create  "smart”  response  surfaces.  In  contrast  to  the  process  independent  polynomial  fit 
of  the  conventional  response  surface  method,  smart  response  surfaces  derive  their 
basic  shape  from  the  process  physics  and  are  then  calibrated  using  designed 
experiments.  This  method  provides  for  a  surface  of  better  representational  accuracy 
using  the  same  or  fewer  experimental  points. 


This  method  has  been  applied  to  the  development  of  a  model  for  the  low  pressure 
chemical  vapor  deposition  (LPCVD)  of  polysilicon,  a  process  used  in  the  manufacture  of 
VLSI  circuits.  A  one-dimensional  finite  difference  model  of  the  LPCVD  process  was 
constructed.  A  Taguchi  orthogonal  array  experiment  was  conducted.  A  confirming 
experiment  performed  at  the  parameter  levels  indicated  by  the  Taguchi  optimization, 
served  to  confirm  the  validity  of  the  experimental  procedure.  The  experimental  results 
will  subsequently  be  used  to  calibrate  the  mechanistic  model. 
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ABSTRACT 


A  meihodology  is  presented  for  the  construction  of  process 
models  by  the  combination  of  physically  based  mechanistic 
modeling  and  statistical  experimental  design  in  order  to  create 
"smart"  response  surfaces.  In  contrast  to  the  process  independent 
polynomial  fit  of  the  conventional  response  surface  method,  smart 
response  surfaces  derive  their  basic  shape  from  the  process  physics 
and  ate  then  calibrated  using  designed  experiments.  This  method 
provides  for  a  surface  of  belter  representational  accuracy  using  the 
same  or  fewer  experimental  points. 


This  method  has  been  applied  to  (he  development  of  a 
model  for  the  low  pressure  chemical  vapor  deposition  (LPCVD)of 
polysilicon,  a  process  used  in  the  manufacture  of  VLSI  circuits.  A 
one-dimensional  finite  difference  model  of  the  LPCVD  process 
was  constructed  A  Tagitchi  orthogonal  array  experiment  was 
conducted  A  confirming  experiment  performed  at  the  parameter 
levels  indicated  by  the  Taguchi  optimization,  served  to  confirm  the 
validity  of  the  experimental  procedure.  The  experimental  results 
will  subsequently  be  used  to  calibrate  the  mechanistic  model. 
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MOTIVATION  FOR  PROCESS  MODELING 


In  the  most  general  sense,  a  process  model  is  a  body  of 
knowledge  which  provides  predictions  about  the  outputs  from  a 
manufacturing  process  given  information  about  inputs  to  the 
process.  Figure  I  illustrates  a  generic  manufacturing  process  and 
related  processing  equipment  The  inputs  have  been  divided  into 
usefully  distinct  categories;  process  parameters  or  control  factors 
and  noise  factors.  Hie  process  parameters  are  those  parameters 
that  wc  cxcicisc  direct  control  over  Examples  of  process 


parameters  include  temperature,  pressures,  roll  speeds,  gas  flow  % 
rates,  etc.  A  second  set  of  inputs  is  entitled  disturbances  or  nois^VV ■ A 
factors.  Noise  factors  are  inputs  lo  the  process  which  are  subject 
unintended  and  undesired  variation.  Examples  of  noise  factors  * 
include  variations  in  the  properties  of  incoming  raw  material,  and 
in  the  process  parameters  themselves.  The  goal  of  the  process 


in  the  process  parameters  themselves.  The  goal  of  the  process 
model  is  to  provide  information  about  the  output  from  tlie  process 
given  information  about  the  process  parameters  and  the  noise 
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Figure  1.  Representation  of  a  Generic  Manufacturing  Process 


Process  modeling  has  several  critical  roles  to  play  A 
competent  model  of  a  manufacturing  process  is  an  extremely 
powerful  tool  in  the  design  of  new  processing  equipment,  as 
accurate  predictive  capability  will  substantially  reduce  the  number 
of  iterations  necessary  to  achieve  a  satisfactory  design. 


Process  modeling  is  also  essential  to  the  operation  of 
existing  equipment.  Models  can  be  used  to  optimize  the  operation 
of  existing  equipment.  Optimization  might  lie  accomplished  by  , 
selection  of  a  set  of  process  parameters  which  leads  to  the  grcatci***- 
robustness  of  the  process  against  disturbances  or  noise  factors 
Process  models  are  also  extremely  useful  for  on-line  quality 
control.  While  we  may  think  of  process  optimization  as  the 
selection  of  an  operating  point  for  the  process,  on-line  quality 
control  deals  with  process  operation  around  that  operating  point. 
Models  can  be  used  to  effectively  guide  the  process  back  to  target 
values. 


The  construction  and  utilization  of  process  models  is 
especially  critical  in  modern  manufacturing  environments 
Competitive  pressures  dictate  that  processes  must  be  run  near  their 
optimum  conditions.  Computer  integrated  manufacturing  offers  a 
potential  wealth  of  data  from  process  operations  w  hich  can  only  be 
effectively  utilized  in  combination  with  process  models 


APPROACHES  TO  MODEL  CONSTRUCTION 


Models  may  be  constructed  by  three  distinctly  different 
approaches;  experimental,  experiential,  and  analytical. 


Experimentally  based  process  models  are  constructed  by 
performing  deliberate  and  planned  experiments  on  a  pintcx* 
Ihese  experiments  are  most  effectively  performed  using  techniques 
of  statistical  experimental  design.  In  these  methods  such  as  Box 


"factorial  experimental  design  and  response  surfaces"  |l,  2|  and 
Taguchi  "orthogonal  array*  [3,  4,  5,  6],  many  experimental 


parameters  are  varied  simultaneously  in  a  well-defined  plan, 
resulting  in  great  economy  of  experimentation.  These  methods  aic 
reviewed  briefly  in  a  later  portion  of  this  paper. 


Experiential  process  models  may  be  constructed  simply  b> 
operating  the  process  in  production  and  collecting  and  analyzing 
the  "observational"  or  "happenstance"  data  that  results  f 


Analytical  models  are  based  on  the  fundamental  physical 
mechanisms  of  the  problem.  Analytical  models  can  l>c  cither 
closed  form,  or  numerical  methods. 
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Physically  based  analytical  models  offer  the  advantage  of 
having  the  greatest  generality  and  range  of  application,  however, 
.  they  CM  be  extremely  time  consuming  to  develop,  and  are  often  of 
•^•!!fuestion,b,e  *ccur,cy  “d  u*c  because  of  a  lack  of  complete 
.■^.knowledge  about  the  process  physics.  Experimental  and 
expenentiil  models  offer  the  advMtage  of  good  fidelity  within  the 
range  of  variables  tested  but  limited  extrapolation  capacity  beyond 
that  range  and  limited  extension  to  equipment  other  than  that  upon 
which  the  experiments  were  run. 


In  today’s  practice,  the  three  methods  of  model  building, 
experimental,  experiential,  Md  analytical,  are  generally  applied 
inttependently  with  little  interaction  between  the  methodologies. 
This  papers  concerns  the  fusion  of  Malytical  and  experimental 
modeling  in  order  to  gain  the  generality  of  m  Malytical  model  in 
combination  with  the  precision  Md  ease  of  use  of  an  experimental 
model. 


BACKGROUND  -  DESIGNED  EXPERIMENTS 


Experimental  design  is  a  systematic  Md  organized  way  to 
conduct  experiments  in  order  to  extract  the  maximum  information 
from  the  minimum  number  of  experiments.  The  unifying  feature 
of  statistically  designed  experiments  is  that  all  the  parameters  of 
interest  are  varied  simultaneously,  in  contrast  to  the  more 
conventional  one  variable  at  a  time  experimental  technique.  In  this 
manner,  the  total  experimental  range  is  explored  with  a  minimum 
number  of  experiments.  There  are  two  commonly  used 
methodologies  for  experimental  design:  the  Taguchi  orthogonal 
array  method  and  the  Box  response  surface  method. 


In  the  Taguchi  orthogonal  array  method,  scientific  and 
engineering  knowledge  is  used  to  pick  experimental  parameters 
which  are  non-interacting.  Typically,  three  levels  would  be 
assigned  for  each  parameter,  low,  medium  and  high.  The 
parameters  or  "control  factors"  are  then  arranged  in  an  orthogonal 
gc'  '-■'xperimental  array.  Figure  2  shows  a  four  parameter,  three  level, 
a)  •jrthogonal  array  which  defines  nine  experiments.  Also  shown  in 
Figure  2  is  a  graphical  representation  of  the  distribution  of 
experimental  points  for  three  levels  of  three  control  factors  or 
parameters.  The  points  plotted  correspond  to  the  second,  third  and 
fourth  columns  of  the  orthogonal  array  shown  Such  a  plot  is 
useful  in  visualizing  the  distribution  of  experimental  points  in 
space,  but  loses  its  utility  past  three  control  factors.  The  unique 
feature  of  m  orthogonal  array  is  that  for  a  given  level  of  a  given 
parameter,  all  other  levels  of  all  other  parameters  are  explored 
uniformly.  Thus,  for  example,  in  runs  one.  two  and  three  of  the 
orthogonal  array  in  Figure  2,  parameters  two,  three  Md  four  are  all 
routed  through  their  low,  medium  Md  high  values. 
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Figure  2.  A  Taguchi  Orthogonal  Array  Md  3  dimensional  • 
visualization  of  columns  2,  3  and  4. 


Another  uniaue  feature  in  the  Taguchi  method  is  the 
ineous  use  of  "Inner”  and  "outer”  arrays.  The  arrav  shown  in 


simultaneous  use  of  "inner”  and  "outer"  arrays.  The  array  shown  in 
Figure  2  is  m  inner  array,  and  is  used  to  define  the  points  of 
interest  in  experimental  space.  An  outer  array  would  be  used  to 
define  small  regions  around  each  of  the  points  specified  in  the 
inner  array.  The  outer  array  would  specify,  for  example,  variations 
or  tolerances  on  experimental  parameters  or  other  "noise"  factors, 
thereby  defining  additional  experiments  in  the  neighborhood  of 
each  of  the  nine  experiments  specified  by  the  inner  array. 


Perhaps  the  most  distinct  feature  of  the  Taguchi  method  is 
the  interpreution  of  tesulu  by  using  a  "signal  to  noise  ratio".  This 
method  of  interpretation  will  be  discussed  in  greater  detail  later  in 
the  paper. 


The  alternative  method,  developed  by  George  Box  and 
others,  defines  a  aeries  of  experiments  and  summarizes  the  results 
of  those  experiments  in  the  form  of  a  response  surface.  A  response 
surface  is  a  polynomial  fit  (usually  a  quadratic  polynomial)  to  the 
measured  data.  The  concept  of  a  response  surface  and  its 
Malytical  representation  is  shown  in  Figure  3  for  a  function  of  two 
variables,  xi  and  X2.  Since  the  response  surface  includes  the  effect 
of  factor  interactions,  a  larger  number  of  experiments  is  needed  to 
fit  a  response  surface  as  compared  to  the  Taguchi  orthogonal  array 
method.  For  example,  a  second  degree  polynomial  fitted  to  four 
parameters  tt  three  levels,  would  require  s  minimum  of  fifteen 
experiments  as  opposed  to  the  nine  experiments  used  in  the 
Taguchi  method. 
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Figure  3.  An  Illustration  of  a  Box  Response  Surface 


The  response  surface  method  is  powerful  in  its  generality  , 
but  suffers  from  the  fact  that  it  does  not  directly  embody  the 
sensitivity  of  the  output  or  response  to  small  deviations  of  input 
f acton.  Since  the  response  surface  is  a  quadratic  fit  to  three  points, 
one  cannot  expect  that  the  slope  at  the  three  points  is  particularly 
accurate.  Since  it  is  the  local  slopes  that  embody  the  sensitivity 
infoimatioa  which  will  be  necessary  to  design  a  robust  process,  the 

Suadrahc  fit  response  surface  is  not  particularly  useful  for  the 
esign  of  robust  processes. 


design  of  robust  processes. 

BUILDING  "SMART"  RESPONSE  SURFACES 


The  beak  conceptual  framework  of  the  present  work  is  to 
replace  the  polynomial  fit  of  a  response  surface  with  a  shape 
dictated  by  the  physics  of  the  process  under  study.  Since  the 
general  shape  will  be  dictated  tty  the  physics,  a  more  precise  model 
can  be  obtained  with  the  same  number  of  experimental  points,  and 
perhaps  even  with  fewer  experimental  points.  In  addition  to  being 
more  accurate  within  the  experimental  range,  such  a  model  would 
also  be  mote  useful  when  extrapolated  beyontj  the  tested  range. 


I 


The  remainder  of  this  paper  describes  one  method  of 
constructing  such  i  "smart”  response  surface  The  steps  in 
construction  are  as  follows: 

*  identify  and  characterize  parameters  or  factors 

*  perform  designed  experiments 

*  develop  a  simple  analytical  model 

*  calibrate  the  model 

*  use  the  model  for  a  design  and  operation  of  equipment 

BACKGROUND  -  LOW  PRESSURE  CHEMICAL  VAPOR 
DEPOSITION  OF  POLYSILICON 

The  process  modeled  in  this  paper  is  a  low  pressure 
chemical  vapor  deposition  process  used  in  the  semiconductor 
industry  to  fabricate  VLSI  circuits.  Integrated  circuits  are  basically 
fabricated  by  alternately  depositing  and  selectively  removing  layers 
on  a  silicon  wafer.  In  the  dominant  family  of  CMOS  circuits,  an 
individual  transistor  has  three  contacts;  the  source,  drain,  and  gate, 
as  shown  in  Figure  4  [7],  The  most  commonly  used  material  for 
the  gate  electrode,  is  polycrystalline  silicon,  which  is  deposited  on 
the  wafer  approximately  midway  through  the  fabrication  of  a 
CMOS  circuit.  Polysilicon  is  used  as  a  conductor  at  this 
intermediate  fabrication  stage,  because  it  allows  the  wafer  to  be 
exposed  to  subsequent  high  temperature  processing  which  metallic 
contacts  would  preclude. 


Figure  4.  Schematic  Cross  Section  of  a  MOSFET  Transistor. 

Polysilicon  layers  are  most  commonly  deposited  on  wafers 
in  a  batch  operation  in  a  tube  furnace  as  shown  in  Figure  5  [8,  9], 
A  tube  furnace  typically  consists  of  a  quartz  tube  sun-ounded  by 
resisting  heating  coils  which  are  in  turn  surrounded  by  thermal 
insulation.  The  four  to  six  inch  silicon  wafers  are  held  25  at  a  time 
in  "boats".  Four  to  six  boats  are  loaded  into  the  tube  furnace  for 
each  batch.  After  loading,  a  partial  vacuum  is  applied  to  the  tube 
and  a  process  gas  in  introduced  through  small  tubular  "injectors  ’. 
For  low  pressure  chemical  vapor  deposition  (LPCVD)  of 
polysilicon,  the  process  gas  is  S1H4.  At  the  operating  temperature 
of  approximately  625°  C„  SiH4  pyrolytically  decomposes  to  yield 
solid  silicon  which  deposits  on  all  hot  surfaces,  and  gaseous 
hydrogen,  according  to  the  relationship: 

SiH4-623°C>  Si  ♦  2H2 


Figure  S.  A  tube  furnace  used  for  the  deposition  of 
polycrystalline  silicon. 


As  the  reacnon  proceed*,  the  liberated  hydrogen  joins  the 
process  gas  in  the  tube,  thus  diluting  the  mixture  and  reducing  the 
reaction  rate.  It  is  primarily  this  dilution  which  leads  to  non- 
uniformitie*  in  deposition  or  growth  rate  from  wafer  to  wafer  down 
the  length  of  the  tube.  Wafer  to  wafer  non-uniformity  is  fVs 
primary  source  of  concern  in  LPCVD  of  polysilicon  For  reascV** 
explained  in  detail  below,  across  the  wafer  uniformity  tends  to  be  ' 
quite  good  in  this  process.  In  order  to  counteract  the  S1H4 
depletion  and  improve  the  wafer  to  wafer  uniformity  in  the  tube. 
SiH4  is  introduced  or  "injected"  at  three  sites  in  a  typical  tube 
furnace  as  indicated  schematically  in  Figure  5. 

Process  control  is  achieved  by  adjustment  of  the  flowrates 
through  each  of  the  three  injectors,  adjustment  of  the  axial  position 
of  the  injectors,  temperature  adjustment,  and  control  of  tube 
pressure. 

IDENTIFICATION  AND  CATEGORIZATION  OF 
PARAMETERS 

The  primary  control  factors  or  parameters  are  tube  pressure, 
flowrate  through  each  injector,  position  of  each  injector,  and 
temperature  profile  down  the  tube. 

The  possible  list  of  noise  factors  is  quite  long,  but  might 
prominently  include  such  issues  as  the  amount  of  prior  deposition 
in  the  tube,  the  prior  condition  of  the  wafers,  aging  of  the  control 
thermocouples,  spacing  between  the  wafers,  and  location  accuracy 
of  the  wafers  in  the  tube. 

The  output  or  response  factor  chosen  for  the  current  model 
is  the  deposition  thickness  as  a  function  of  position  down  the  tube 
In  practice,  one  is  also  concerned  with  controlling  the  grain  size  of 
the  deposition,  and  the  thickness  distribution  across  each  wafer. 

As  noted  earlier,  the  thickness  across  a  wafer  tends  to  be  quite 
uniform.  As  the  grain  size  is  controlled  almost  exclusively  by 
deposition  temperature,  the  desired  grain  size  fixes  the  temperature 
of  deposition  typically  at  625“  G 

For  the  purposes  of  our  modeling,  we  have  chosen  the  four 
control  factors:  tube  pressure,  the  flowrate  through  the  load  end 
injector,  the  flowrate  through  the  center  injector,  and  the  position 
of  the  source  injector.  The  load  end  injector  is  the  injector  on  the 
left  side  of  Figure  5.  The  source  end  injector  is  the  injector  whose 
opening  within  the  tube  furnace  is  furthest  to  the  right  in  Figure  5 
The  output  or  response  factor  in  our  experiment  will  be  the  profile 
of  thickness  down  the  length  of  the  tube  furnace. 

EXPERIMENTAL  DESIGN.  RESULTS.  AND 
INTERPRETATION 

The  experimental  design  used  was  a  Taguchi  orthogonal 
array  using  four  parameters  at  three  levels.  This  array,  shown  m 
general  form  in  Figure  2,  is  again  shown  in  Figure  6  complete  with 
parameter  assignments  and  level  selections.  Three  of  the  four 
pvameters  (all  but  pressure)  are  indicated  in  dimensionless  form 
The  nine  experiments  were  conducted  with  a  wafer  load  of  1 50  si  x- 
inch  wafers.  Thirteen  test  wafers  were  distributed  within  the  150 
wafer  load.  A  baseline  experiment  was  repeated  five  times  in  order 
10  gain  some  information  about  run  to  run  variability  of  the 
process. 

Figure  7  shows  plots  of  growth  rate  (avenged  over  the  test 
wafer)  against  position  in  the  tube  furnace  for  runs  one  and  nine  of 
the  orthogonal  array.  These  plots  are  indicative  of  the  range  of 
results  obtained.  The  mean  values  and  standard  deviations  of  each 
of  the  thirteen  test  wafer  positions  are  shown.  The  mean  value  has 
been  obtained  from  a  single  run  of  each  experiment.  The  standard 
deviation  was  obtained  by  normalizing  the  standard  deviation  for 
each  wafer  position  in  the  baseline  replicate  tuns  and  applying  this 
normalized  standard  deviation  to  the  nine  Taguchi  arrj^. 
experiments.  The  bar  chans  at  the  bottom  of  the  plots  in  Figur..>\.^ 
schematically  indicates  the  gas  flow  at  each  injector  site.  As  m.’.-' 


be  seen,  there  is  a  general  correlation  between  local  injector 
volume  and  growth  rate.  The  system  thus  seems  to  have  no 
conspicuous  pathologies. 
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Figure  6.  Taguchi  Orthogonal  AtTay  Experimental  design  used. 

Qload  +  Qcenter  +  Qsource  *  150  stnd  cm^/min. 

Process  temperature  *  625°. 

The  results  of  the  designed  experiments  are  to  be  used  in 
two  ways.  Later  they  will  be  used  to  calibrate  the  smart  response 
surface  model.  First,  however,  they  will  be  used  to  predict  an 
"optimized"  point  of  operation  of  the  equipment.  This  predicted 
optimum  will  then  be  run  in  a  confirming  experiment.  If  the 
resulting  improvement  is  close  to  that  predicted,  we  may  be 
confident  that  the  experimental  parameters  and  output  variables 
have  been  properly  identified  and  that  the  experiments  were 
if  .■performed  well. 

*  -  The  interpretation  of  the  orthogonal  anay  results  is 

indicated  schematically  in  Figure  8.  and  begins  by  calculation  of  a 
signal  to  noise  ratio  for  each  of  the  nine  experimental  runs.  This 
signal  to  noise  ratio  characterizes  the  deviation  of  each  of  the 
profiles  from  a  flat  and  uniform  profile.  Next,  average  signal  to 
noise  ratios  are  calculated  for  each  level  of  each  parameter,  and  are 
plotted  on  the  marginal  graphs  of  Figure  9.  These  graphs  may  now 
be  used  to  select  the  best  combination  of  parameter  levels  for  an 
optimized  process.  This  optimum  occurs  at  middle  values  for  load 
injector,  flow  source  injector,  low  values  for  flow  source  injector 
position  and  pressure. 

This  combination  of  parameters  may  now  be  run  in  a 
confirming  experiment,  the  results  of  which  are  shown  in  Figure 
10.  As  may  be  seen,  the  confirming  experiment  is  a  substantial 
improvement,  thus  validating  the  experimental  procedure. 


DEVELOPMENT  OF  THE  ANALYTICAL  MODEL 

The  development  of  the  analytical  model  begins  by 
performing  order  of  magnitude  calculations  in  order  to  identify 
important  physical  mechanisms.  First,  mass  transport  calculations 
indicate  that  the  the  deposition  rate  is  limited  by  the  surface 
reaction  rate,  and  not  by  transport  to  the  surface.  The  flow  in  the 
annular  region  between  wafers  and  tube  is  characterized  by  a 
Reynold's  number  of  approximately  one,  indicating  that  the  flow  is 
laminar.  The  Pedet  number  in  the  annular  region  is  approximate  y 
one.  indicating  that  convective  and  diffusive  fluxes  are  roughly 
comparable,  and  that  both  must  be  considered  in  the  solution  to  the 
<Tr’1problcm.  A  Poiselle  flow  calculation  in  the  annular  space  indicates 
the  pressure  drop  down  the  length  of  the  tube  is  less  than  one 
"'‘‘'"'‘percent  of  the  actual  pressure  in  the  tube,  and  hence  is  negligible. 
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Figure  7.  Measured  growth  rate  P‘^»«"ns!  Wu,t  f°r 
8  wo  of  the  nine  experiments  defined  in  Figure  6 

In  the  region  between  wafers,  the  Peclet  number  is  much  less  than 

one.  indicating  that  diffusive  nu*e*^in*^*enxdit'velocilv 
convective  flow  in  this  region  can  be  ignored.  The  exi  i  velocity 

from  the  small  diameter  injector  tubes  is  °"J* 

1  and  hence,  the  flow  geometry  and  incorpora  ion  into  ihe  tube 
annular  space  is  quite  complex.  The  gas  depletion,  that  is  e 

percentage  of  SiHg  that  is  reacted  to  form  silicon,  is  between  20*1 
and  50%,  and  is  thus,  an  important  component  of  the  problem 
Indeed,  it  is  this  depletion  that  is  the  primary  source  of  the 
variation  in  deposition  rate.  The  surface  reaction  rate  is  a  function 
of  the  local  partial  pressures  of  SiHa  and  hydrogen  and  the  local 
temperature. 

Consistent  with  these  order  of  magnitude  calculation,  our 
model  assumes  a  one  dimensional  finite  difference  formulation 
with  no  radial  non-uniformities.  The  flow  has  been  assumed  to  be 
inviscid,  and  to  be  of  a  piug  flow  nature.  We  are  concerned  with 
the  coupled  problem  of  convection  and  diffusion  in  the  annular 
space.  We  will  consider  the  thermal  problem  to  be  separated  from 
the  flow  problem  and  will  specify  the  wafer  temperature  as  a 
function  of  position  down  the  tube.  We  will  also  specify  the 
annular  flow  area  as  a  function  of  position  down  the  tube,  we  ill 
flwtoi  the  injector  flow  as  mixing  into  the  annular  space  in  the  tube 
over  a  tube  length  that  is  proportional  to  the  square  of  the  flowrate 
(and  velocity)  tit  the  tube. 
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Figure  8.  The  experimental  results  are  interpreted  by  calculating  a 
signal  to  noise  ratio  for  each  of  the  nine 
experiments  of  the  orthogonal  array. 
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Figure  9.  The  "marginal  graphs"  which  plot  avenge 
signal/noise  ratio  for  each  level  of  each  parameter. 


WAFER  POSITION 


Figure  10.  Confirming  experiment  at  "optimized"  control  factors 

The  model  basically  consists  of  mass  conservation 
conditions  for  each  of  the  two  species,  Sitt*  and  hydrogen  taken 
individually  and  together  as  summarized  in  the  equations  below. 

Mass  Conservation  of  Silane: 


±(AcD£- 

dz  v  dz 


AcVx,  «  R(x,)-F 


Mass  Conservation  of  both  species: 


—  (cAV)  »  R (x,)  +  F 
dz 

where: 

A  is  the  cross  sectional  area  of  the  flow  [cm2]. 

C  is  the  total  concentration  of  gaseous  species  [moles/cm3], 
z  is  the  distance  down  the  tube  [cm], 

D  is  the  diffusion  of  silane  [cm2/sec], 

x$  is  the  molar  fraction  of  silane  [moles/mole], 

V  is  the  molar  average  velocity  [cm/sec], 

R  is  the  surface  reaction  rate  and  is  a  function  of  xs  and 
temperature,  [mole/sec  -  cm], 

F  represents  the  injector  inlet  flow  [mole/sec  •  cm] 

MODEL  CALIBRATION 

The  final  step  in  construction  of  a  smart  response  surface  is 
the  use  of  the  results  of  the  experimental  design  to  perform  a  least 
squares  fit  to  calibrate  the  model  parameters.  The  adjustable 
parameters  in  our  model  include  four  parameters  that  govern  the 
deposition  rate  as  a  function  of  partial  pressures  and  temperature, 
and  two  parameters  that  specify  the  nature  of  the  incorporation  of 
injector  flow  into  the  annular  space. 

At  the  current  writing,  the  model  has  been  developed  an^xi 
found  to  converge  to  proper  solutions  rapidly.  In  the  near  futurvvyy 
the  model  will  be  calibrated  using  the  experimental  data.  '  ' 
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CONCLUSIONS  AND  FUTURE  WORK 


We  have  illustrated  a  methodology  for  the  construction  of 
process  models  using  the  concept  of  the  "smart'  response  surface. 
A  smart  response  surface  is  one  in  which  the  general  shape  of  the 
response  surface  is  dictated  by  the  process  physics  ant)  the  shape  is 
then  calibrated  using  statistically  designed  experiments.  This 
approach  contrasts  with  the  conventional  response  surface  methods 
in  which  a  general  polynomial  is  fit  to  the  data  derived  by  designed 
experiments.  The  motivation  behind  constructing  a  smart  response 
surface  is  that  the  resulting  model  will  have  greater  fidelity  in 
representing  the  process  both  within  the  range  of  experimentation 
and  when  extrapolated  beyond  the  experimental  range. 

The  smart  response  surface  approach  has  been  discussed 
within  the  context  of  building  a  model  for  the  low  pressure 
chemical  vapor  deposition  of  polysilicon  as  performed  in  the 
integrated  circuit  industry.  In  the  work  discussed,  a  finite 
rlilfercnce  model  of  the  process  was  built,  and  a  Taguchi 
orthogonal  array  experimental  design  was  performed.  The  results 
from  the  designed  experiment  were  used  to  optimize  the  process 
and  a  confirming  experiment  at  the  predicted  optimum  conditions 
demonstrated  the  validity  of  the  experimental  program  In  the 
future,  the  designed  experiments  will  be  used  to  calibrate  the 
numerical  model 

Future  work  will  focus  on  ways  to  combine  designed 
experiments  with  fragmentary  mechanistic  modeling.  This  will 
permit  the  use  of  the  smart  response  surface  technique  without  the 
requirement  that  a  complete  physical  model  be  available. 
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